Spin Wave Propagation in the Domain State of a Random Field Magnet 
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Inelastic neutron scattering with high wave-vector resolution has characterized the propagation 
of transverse spin wave modes near the antiferromagnetic zone center in the metastable domain 
state of a random field Ising magnet. A well-defined, long wavelength excitation is observed despite 
the absence of long-range magnetic order. Direct comparisons with the spin wave dispersion in the 
long-range ordered antiferromagnetic state reveal no measurable effects from the domain structure. 
This result recalls analogous behavior in thermally disordered anisotropic spin chains but contrasts 
sharply with that of the phonon modes in relaxor ferroelectrics. 

PACS numbers: 75.30.Ds, 75.50.Lk 
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Disorder in a condensed matter system can strongly 
influence mode propagation through the medium, pro- 
viding a unique perspective on its microscopic behavior. 
To gain a more comprehensive understanding of such ef- 
fects from disorder, we have performed a neutron scat- 
tering study of the long wavelength spin waves of the 
frozen domain state of a random field magnet. Despite 
the wealth of attention directed at the random field prob- 
lem , little attention has been directed toward the spin 
wave dispersion in these systems. Our study has been 
motivated in part by experiments on the phonons in sev- 
eral relaxor ferroelectrics including Pb(Zn 1 / 3 Nb 2 /3)03 
and Pb(Mg!/3Nb2/3)03 that have demonstrated inter- 
esting anomalies associated with the disorder in these 
systems H, H| - Specifically, at wave vectors near 0.2 
A -1 , the dispersion of the transverse optic branch ap- 
pears to drop precipitously into the acoustic branch. This 
phenomenon has been identified with a sharply wave- 
vector dependent overdamping associated with the pres- 
ence of polar nanoregions in the relaxors . To test the 
generality of such behavior for systems with disorder, we 
have performed a search for analogous effects in the spin 
waves of the domain state of the random field magnet. In 
marked contrast to the phonons in the relaxors, we ob- 
serve well-defined long wavelength excitations transverse 
to the ordering axis, despite the absence of long range 
order, and find no measurable effect on the spin waves 
from the domain structure. We identify the persistence 
of these modes as an apparently generic feature of disor- 
dered magnetic systems with anisotropy, as first realized 
in thermally disordered spin chains 
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The experiments were performed on a diluted antifer- 
romagnet, Mn x Zni_ x F2 with x = 0.5, cooled in a mag- 
netic field parallel to the Ising axis. The physics of the 
diluted Ising antiferromagnet in a uniform field can be 
mapped directly on to that of a random field Ising mag- 
net [9(, providing a method for experimental realization 



of the random field system. MnFa has a tetragonal rutile- 
type structure with lattice constants a = 4.87 A and 
c = 3.31 A. The spin interaction is primarily Heisenberg- 
like, and a weak dipolar interaction between the Mn mo- 
ments accounts for the small Ising anisotropy that aligns 
the Mn spins along the c axis (i.e., the tetragonal axis). 
Diluted MnF2 was chosen for this study because, as a 
weakly Ising antiferromagnet, its small spin gap makes 
the physics of interest easily accessible. We stress that, 
while the interesting phonon behavior in relaxors moti- 
vates our spin wave study on this system, the diluted 
antiferromagnet in field should not be considered a pre- 
cise magnetic analog of disordered ferroelectrics. First, 
the antiferromagnet does not share the property of a con- 
served order parameter. In addition, the nature of the 
polar nanoregions in the relaxors as the source of ran- 
dom fields has recently come into question [ljj. Never- 
the-less, in pursuing the general problem of the effects of 
disorder on the spin waves in a random field magnet, we 
view diluted MnF 2 as the best candidate for experiments. 

When cooled from the paramagnetic phase in the field, 
the diluted antiferromagnet forms a metastable state of 
frozen domains. Using high resolution neutron scatter- 
ing, we have characterized this domain structure and the 
corresponding spin wave dispersion at wavelengths larger 
than the typical domain size. In addition, by cooling the 
antiferromagnet in zero field and then applying the field, 
we have produced a long-range ordered antiferromagnetic 
state in the system. Comparing the results obtained after 
cooling in field with the results of identical measurements 
taken after zero field cooling, we are able to determine 
directly the influence of the domain structure on the spin 
waves. 

An important criterion for this experiment is the abil- 
ity to create domains under field cooling that are suffi- 
ciently small to resolve the spin wave dispersion at the 
appropriate wave vectors. Smaller domains are achieved 
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with higher dilution [Til IT^ . larger applied fields [T^ . 
and faster cooling rates. With too high dilution, how- 
ever, the spin waves (in zero field) eventually become 
poorly defined 01 ■ I n addition, the presence of a spin- 
flop transition in Mn a; Zni_ a ;F2, with a spin-flop field that 
depends on x 0] , limits the maximum field that one can 
apply while remainin g in the random field Ising regime. 
Previous studies [TH TT2I Il3| indicate that x — 0.5 is 
the optimal compromise among these different concerns. 
The experiment was performed at the NIST Center for 
Neutron Research on the SPINS triple-axis spectrometer, 
with tight collimation (30'-10'-S-10'-20') to achieve high 
wave-vector resolution and a final neutron energy fixed 
at 3.5 meV. A single crystal of Mno.5Zno.5F2, approx- 
imately 50 grams in size, was oriented with the c-axis 
(the Ising axis) parallel to the external field and perpen- 
dicular to the scattering plane. The crystal was strongly 
coupled to the bath of a helium flow cryostat, permitting 
rapid cooling (> 20 K/minute) from the high tempera- 
ture paramagnetic phase to 4.2 K, well below the Neel 
temperature in zero field, T^{H = 0) = 21 K. 

Figure 1 displays scattering profiles from elastic scans 
through the antiferromagnetic zone center (1, 0, 0) along 
the transverse direction (l,gx,0) in an external field of 
H = 2.3 T. The open circles represent measurements 
taken after cooling the sample in zero field, so that the 
random field strength is zero through the Ising transition, 
and then raising the field (ZFC). The scattering follows 
a resolution-limited lineshape consistent with long-range 
antiferromagnetic order. The solid circles show the re- 
sults from cooling in field (FC) and display a significantly 
broadened lineshape corresponding to the metastable do- 
main state induced by the random fields. The scattering 
in the field-cooled measurement also shows a large overall 
enhancement in intensity, which we associate with the re- 
lief of extinction due to the domain structure. We model 
the short-ranged correlations of this field-cooled domain 
state with a Lorentzian-squared form, 



where q is the wave vector measured with respect to the 
antiferromagnetic zone center. The solid line in Fig. 1 is 
the results of a fit with Eq. (1) convolved with the in- 
strumental resolution function. This static neutron cross 
section is proportional to the Fourier transform of the 
real-space two spin correlation function. The transform 
of Eq. (1) corresponds to an exponential decay in the 
correlations between Ising spins, 

(5,(r)5,(0))~e-^. (2) 

The Lorentzian-squared form for the elastic cross sec- 
tion has been shown in previous studies of field- 
cooled Mn a; Zni_ a ;F2 [Til IT^ and other random 
field Ising magnets, including Co x Zni_ a; F2 and 
Rb2Co a: Mgi_ 2 ;F4 [15|. to describe accurately the static 
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FIG. 1: Neutron scattering intensity at zero energy transfer 
scanning through the antiferromagnetic zone center (1, 0, 0) 
along the direction (0, q K , 0) at T = 4.2 K and H = 2.3 T. 
The open circles are the results after zero field cooling and 
have a resolution-limited lineshape. The solid circles are the 
results after field cooling and have a lineshape reflecting the 
short-range ordered domain state. The field cooled intensities 
have been divided by 3.5 for comparison. The solid line is a 
fit to Eq. (1) convolved with the instrumental resolution. 



short-range order in the field-cooled state. In addi- 
tion, theoretical arguments have indicated that this form 
should be expected on general grounds 0, 0, 0, 0] . 
Thus, field-cooled state can be considered a collection 
of metastable domains whose size distribution produces 
the real-space correlation function given by Eq. (2). The 
correlation length extracted from the fit with Eq. (1), 
£ = 137 ± 4 A, sets the characteristic length scale in 
the a-b plane for the domains. Assuming that the do- 
main structure is isotropic in units of the lattice spacings, 
as suggested by the crystal structure, the smaller lattice 
constant along the c-axis leads to a characteristic scale 
of 93 ± 3 A in that direction. 

To compare the spin wave behaviors in the long- 
range ordered and domain states, we have measured 
the scattering intensity as a function of energy trans- 
fer, to, at q = through the dispersion, as shown in 
Fig. 2. Because of the finite experimental resolution, 
this measurement effectively integrates over the approx- 
imate wavevector ranges — 0.004A -1 < qjj < 0.004A -1 , 
-0.003A- 1 < q K < 0.003A" 1 within the scattering 
plane and — 0.085A -1 < q^ < 0.085A -1 perpendicular 
to it [20L l2l| . Based on the phonon behavior observed in 
the relaxor ferroelectrics, we expect any change caused 
by the domain structure to affect the spin waves at wave 
vectors less than q ss 2-7r/£ = 0.046 A -1 along H and 
K and q w 2?r/£ = 0.068 A" 1 along L. Thus, a re- 
distribution of spectral weight to lower energy, like the 
dramatic overdamping seen in the relaxors, should be 
clearly visible in the scattering intensity in Fig. 2. As 
the figure illustrates, the results for the long-range or- 
dered antiferromagnet and domain state are essentially 
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FIG. 2: Neutron scattering intensity as a function of energy 
transfer, u, at the antiferromagnetic zone center (1, 0, 0) af- 
ter field cooling (solid circles) and zero field cooling (open 
circles). The essentially identical profiles through the spin 
wave dispersion at lu ~ 0.52 meV demonstrates that the long 
wavelength spin waves are insensitive to the domain struc- 
ture. The solid line is the calculated result for the spin wave 
peak plus background based on the spectrometer resolution 
function and the known dispersion relation in Mno.5Zno.5F2 
[13], assuming the that intrinsic width of the dispersion is 
zero. The enhanced intensity after field cooling at small ui, 
shown in the inset, reflects the enhanced Bragg scattering due 
to relief of extinction. 



identical through the dispersion at lu = 0.5 meV, indi- 
cating that the long wavelength spin waves are insensi- 
tive to the domain structure. The scattering at smaller 
energy transfers, shown in the inset to Fig. 2, does dis- 
play enhanced intensity in the domain state. However, 
we associate this enhancement with the stronger Bragg 
scattering observed for the domain state in Fig. 1, ex- 
tending away from lu = due to the finite experimental 
resolution, and not to any change in the spin wave behav- 
ior. In particular, any redistribution of spectral weight 
should be reflected in both enhanced scattering at small 
lu and reduced scattering at larger lu, which is excluded 
by the results in Fig. 2. This dramatic contrast with 
the phonons in relaxors lends additional support to re- 
cent evidence questioni ng t he accuracy of a random field 
picture for the relaxors 10]. 

The solid line in Fig. 2 is the calculated result for the 
spin wave peak plus background based on the spectrom- 
eter resolution function and the previously measured dis- 
persion relation in Mno.5Zno.5F2 [l3T ]. assuming the that 
intrinsic width of the dispersion is zero. The only free 
parameters entering the calculation are the overall am- 
plitude of the peak and the form of the sloping back- 
ground, Alu~ 2 + B. The strongly asymmetric lineshape 
for the calculated result derives from the non-zero slope 
of the dispersion relation away from the zone center and 
the relatively coarse wave-vector resolution perpendicu- 



lar to the scattering plane |21|. Any small discrepan- 
cies between the calculation and measurements likely re- 
sult from uncertainties in this contribution to the reso- 
lution and in the dispersion relation, implying that the 
long wavelength spin waves in both field cooled and Neel 
states are long-lived. Furthermore, direct comparisons 
of the amplitudes and widths of the measured FC and 
ZFC peaks indicate that they are the same to within 5%, 
setting a small upper bound on any possible broadening 
specific to the field cooled domain state. 

This clear insensitivity of the long wavelength spin 
waves to the static disorder of the FC domain structure 
raises two interesting issues that merit further study. The 
first concerns the effect on short wavelength spin waves. 
A comparison between the spin wave dispersions in the 
FC and ZFC states at wave vectors close the antifer- 
romagnetic zone edge would address this question. As 
shown in Ref . |13| , the spin wave dispersion at wave vec- 
tors away from the zone center in a diluted Ising anti- 
ferromagnet becomes increasingly diffuse due to the ge- 
ometric disorder of the dilution. While this effect could 
severely complicate comparisons between FC and ZFC 
behaviors, it also indicates that short wave length spin 
waves are strongly affected by disorder. The second issue 
for further study concerns the temperature dependence 
of the spin wave dispersion, particularly in the transi- 
tion region. The interplay of the critical behavior and 
the disorder-induced slow dynamics in the random field 
Ising magnet has presented serious challe nges to a com- 
prehensive description of the transition [231,13113 ■ A 
characterization of the evolution with temperature of the 
spin wave dispersion near the zone center on approach to 
the transition would provide a new perspective on this 
issue. 

The energy gap of approximately 0.52 meV for the 
excitation in Fig. 2 matches that expected from the 
Heisenberg exchange field and dipolar anisotropy of 
Mno.5Zno.5F2, approximately 0.27 meV (l3|. and the 
Zeeman energy from the precession around the applied 
field, 0.25 meV for a 2.3 T field. This excitation, the 
gapped mode of the Ising state, corresponds to long wave- 
length oscillations in the spin components transverse to 
the easy axis. Because of the experimental resolution, 
the peak in Fig. 2 represents contributions from oscil- 
lations over a range of wavelengths, most of which far 
exceed the correlation length for magnetic order in the 
domain state. The existence of these modes, despite 
the disorder, bears a striking similarity to the behav- 
ior in thermally disordered spin chains with easy plane 
anisotropy. As first predicted by Villain well-defined 
long wavelength transverse modes have been observed 
in a number of one-dimensional anisotropic spin systems 
including CsNiF 3 and (CD 3 ) 4 NMnCl 3 @. The mea- 
surement of such excitations in the domain state of a 
three-dimensional random field Ising magnet indicates 
that these long-lived transverse modes occur more gen- 
erally than previously realized, irrespective of the na- 
ture of the anisotropy, the system's spatial dimension, or 
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the type of disorder. Indeed, studies of spin glasses ap- 
pear to support this picture that these modes are generic 
to anisotropic disordered magnets. While spin waves in 
canonical Heisenberg spin glasses [2(| are strongly over- 
damped, well-defined long wavelen gth excitations have 
been observed in Ising spin glasses |2jj. Theoretical in- 
sight into the apparently universal nature of these trans- 
verse modes in disordered anisotropic magnets would be 
invaluable. 
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